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1.  INTRODUCTION 


1.1.  The  concepts  of  moment  of  inertia  and  shear  area  are 
fundamental  to  any  discussion  of  the  strength  of  materials  or  to 
the  load  bearing  capability  of  a structure.  In  particular,  the 
design  of  ships  requires  a knowledge  of  these  quantities  in  order 
to  compute  the  longitudinal  strength.  This  is  essential  to  deter- 
mining whether  the  girder  is  sufficiently  strong  to  carry  the 
imposed  weight  of  the  ship,  its  stores  and  on-board  equipment. 
Standard  procedures  for  computing  these  quantities  have  been 
developed.  These  are  described  in  reference  1.  The  purpose  of 
this  NSWC/WOL/TR  is  to  describe  the  general  principles  involved 
and  how  they  have  been  implemented  by  computerization. 

2 .  THEORY 


2.1.  In  longitudinal  strength  calculations  of  a ship,  simple 
beam  theory  is  assumed  to  hold,  reference  2.  That  is  to  say,  a ship, 
for  purposes  of  strength  considerations,  may  be  likened  to  a hollow, 
nearly  rectangular,  girder.  Therefore,  the  equations  derivable  from 
simple  beam  theory  may  be  used  to  make  calculations  for  the  strength 
of  such  a beam  of  known  cross-section,  reference  3.  Because  the 
girder  may  be  considered  to  be  a beam  with  a known  cross-sectional 
area  A,  those  simple  relations  which  may  be  found  in  any  text  book 

on  engineering  mechanics  may  be  applied  to  it;  for  example,  see 
reference  4.  However,  for  convenience,  the  relevant  relations  are 
rederived  in  the  sections  below. 

2.2.  Consider  a beam  with  the  arbitrary  cross-section  shown  in 
Figure  2.1.  The  total  area  A can  be  found  by  integration  as  follows: 


The  distance  y to  the  centroid  C of  the  area  is  obtained  from 


where  it  is  understood  by  the  symbol  A beneath  the  integral  signs 
that  the  integrations  must  be  carried  out  over  the  entire  area  A. 


1.  DDS  9290-8.  "Longitudinal  Strength  Calculations,"  Design  Data 
Sheet,  Department  of  the  Navy,  Bureau  of  Ships,  27  Dec  1950, 
Unclassified. 

2.  'Principles  of  Naval  Architecture,"  revised  edition  published  by 
the  Society  of  Naval  Architects  and  Marine  Engineers,  1967. 

3.  T.  C.  Gillmer,  "Modern  Ship  Design,"  United  States  Naval  Institute, 
Annapolis,  Maryland,  Copyright  1970. 

4.  S.  P.  Timoshenko  and  J.  M.  Gere,  "Mechanics  of  Materials,"  D.  Van 
Nostrand  Company,  Copyright  1972. 
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The  numerator  of  Equation  2.2  is  known  as  the  first  moment  of  the 
area  with  respect  to  the  x -axis  and  will  be  denoted  by  the  symbol 
Q^.  Thus,  ^ 

“x  'X  <2-3) 

2.3.  Before  continuing,  it  is  necessary  to  introduce  the  concept 
of  neutral  axis.  If  a beam  is  bent,  some  of  the  fibers  in  the  beam 
will  be  in  compression  while  others  will  be  in  tension.  However, 
there  will  exist  some  surface  within  the  beam  such  that  all  fibers 
lying  in  this  surface  are  neither  in  compression  nor  in  tension. 

This  surface  is  referred  to  as  the  neutral  surface  and  its  inter- 
section with  any  cross-sectional  plane  of  the  beam  is  called  the 
neutral  axis  of  the  cross-section.  The  moment  of  inertia  of  a cross- 
section  of  a girder  used  in  longitudinal  strength  calculations  is 
the  moment  of  inertia  with  respect  to  the  neutral  axis  of  the  cross- 
section.  The  moment  of  inertia  of  an  area  with  respect  to  an  axis 
in  the  plane  of  the  area  is  related  to  the  moment  of  inertia  with 
respect  to  a parallel  centroidal  axis  by  the  Parallel  Axis  Theorem. 

To  obtain  this  theorem,  consider  the  area  shown  in  Figure  2.1. 

Assume  that  the  point  C is  the  centroid  of  the  area  and  that  the  x 
and  y axes  are  through  the  centroid.  The  Xj^  and  y^^  axes  are  parallel 
to  the  X and  y axes  and  have  their  origin  at  any  point  0.  Then,  by 
definition,  the  moment  of  inertia  of  the  area  with  respect  to  the  x^^ 
axis  is 


I 


X 


(y 


+ y)^  da 


or. 


(2.4) 


I 


X 


da  + 


yda  + f 

A 


(2.5) 


The  first  integral  on  the  right  hand  side  of  Equation  2.5  is  the 
moment  of  inertia  about  the  centroidal  x-axis,  the  second  integral 

vanishes  because  the  x-axis  passes  through  the  centroid; 
and  the  third  integral  is  the  area  A of  the  figure.  Therefore, 
Equation  2.5  may  be  reduced  to  the  form 


Ix  = ^c  y^  A.  (2.6) 

It  will  not  be  proved  here  but  it  may  be  shown,  reference  4,  that  the 
neutral  axis  passes  through  the  centroid,  hence,  if  Ij^^  represents 
the  moment  of  inertia  with  respect  to  the  neutral  axis,  we  have  from 
Equation  2.6, 
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= I -y^  A. 

NA  X 


:2.7) 


This  equation  may  be  used  to  find  the  moment  of  inertia  with  respect 
to  the  neutral  axis  given  the  moment  of  inertia  with  respect  to 
the  arbitrary  x -axis  where 


I = 

X 


da . 


(2.8) 


2.4.  The  procedure  to  be  used  is  the  following.  Choose  any 
arbitrary  axis  parallel  to  the  baseline  of  the  ship  girder  cross- 
section  of  interest  (see  Figure  3.1).  Let  this  axis  be  the  x^^-axis. 
Indeed,  the  baseline  itself  may  be  chosen  as  this  arbitrary  axis. 
Now,  compute  the  cross-sectional  area  A from  Equation  2.1,  the  first 
moment  of  the  area  about  this  axis  from  Equation  2.3  and  finally 
the  moment  of  inertia  from  Equation  2.8.  The  distance  y between 
this  arbitrary  axis  and  the  neutral  axis  may  be  found  from  Equation 
2.2  which,  for  simplicity,  may  be  rev/ritten  as 


y = Q^/A. 


(2.9) 


Finally,  use  Equation  2.7  to  obtain  the  desired  moment  about  the 

neutral  axis. 

2.5.  The  calculation  of  the  area  of  a cross-section  effective 
in  shear,  Ag , is  a difficult  problem  where  the  cross-section  does  not 
have  a geometrically  simple  shape.  Such  is  the  case  for  ship  girder 
cross-sections.  Several  formulae  for  computing  the  shear  area  of 
ship  cross-sections  have  been  proposed  in  the  literature.  Hides, 
reference  5,  discusses  these.  For  our  purposes,  the  following 
formula,  derived  from  considerations  of  simple  beam  theory,  is  used 
to  obtain  the  shear  area,  namely 


A = 
s 


"na  ^ 

M 


(2.10) 


The  quantity  is  obtained  from  Equation  2.7,  M is  the  first 

moment  of  the  area  above  the  neutral  axis,  and  t is  the  sum  of  the 
thicknesses  of  all  continuous  side  shell-  plating  intersected  by 
the  neutral  axis.  By  continuous,  it  is  meant  that  the  plating 


5.  A.  N.  Hicks,  "The  Theory  of  Explosion  Induced  Ship  Whipping 
Motions."  NCRE  Report  R579,  (1972). 


5 


NSWC/WOL/TR  77-26 


extends  longitudinally  at  least  1/20  of  the  water  line  length  of 
the  ship.  As  an  illustration,  if  we  refer  to  Figure  3.1,  we  see 
that  the  neutral  axis  of  the  cross-section  of  the  girder  intersects 
the  E-strake  which  we  know  from  Equation  3.1  to  be  0.5  inches  thick. 
But  the  neutral  axis  intersects  the  E-strake  of  both  sides  of  the 
girder,  hence  t is  the  sum  of  the  thicknesses  of  both  E-strakes 
or  1 inch.  The  quantity  M may  be  calculated  from 


where  the  integral  extends  only  over  the  area  A'  of  the  cross- 
section  above  the  neutral  axis. 


3.  READING  SKIPS'  PLANS 


3.1.  In  order  to  carry  out  the  calculations,  it  is  necessary 
to  have  suitable  ship's  plans.  Such  plans  are  provided  by  inertia 
sections.  These  are  cross-sections  of  the  girder  which  show  only 
those  structural  e’  ents  which  effectively  contribute  to  the 

These  are  typically  provided  for  up  to  15 
D.  Figure  3.1  is  a typical  inertia  section 
ship.  The  reading  of  inertia  sections 
cf  terms  used  by  naval  architects.  Some  of 
these  are  discussed  below. 


longitudinal  str 
stations  along 
taken  at  stat' 
requires  a kn 
the  more  commit 


3.2.  The  following  notations  are  commonly  encountered: 

L - refers  to  an  inner-bottom  longitudinal.  These  are 
normally  sequentially  numbered  from  the  keel, 
e. g. , L-1 , L-2 , etc. 

CVK  - Center  vertical  keel.  A major  strength  element. 

P and/or  S - port  and/or  starboard.  Used  to  mark  those 

items  found  only  on  the  port  and/or  starboard 
side  of  the  ship  cross-section.  An  area  of 
plating  enclosed  by  a triangle  marked  P 
and/or  S means  that  it  is  not  effective  for 
strength  purposes  and  should  be  omitted 
from  strength  calculations. 

A,B,C,  etc.  - these  are  used  to  label  the  various  plating 
elements  (commonly  referred  to  as  strakes) 
that  make  up  the  shell  (outer  portion  of 
hull)  and  the  inner-bottom.  These  are 
labeled  sequentially  from  the  keel.  The 
strake  labeled  "G"  in  Figure  3.1  is  often 
called  the  sheer  strake  and  the  deck  (main) 
plate  connected  to  it  is  called  the  stringer 
strake . 
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FK  - Flat  keel.  A flat  plate  welded  to  the  CVK  e.xtending 
along  the  length  of  the  ship. 

I:  - Center-line.  The  a.xis  of  symmetry  of  the  cross-section. 

^ - denotes  a welded  seam.  These  define  the  extent  of 
various  strakes. 

Plating  thicknesses  are  given  in  pounds,  denoted  # per  square  foot 
of  plating.  Using  the  density  of  steel,  the  thickness  of  a given 
plate  whose  weight  is  W #/ft  sq  is. 


thickness  (in)  = W/40.8.  (3.1) 


For  example,  the  A strake  of  the  inner-bottom  which  entends  from 
the  top  of  the  CVK  to  the  first  welded  seam  is  labeled  "A"  - 20.4  #. 
Its  thickness  is  therefore  1/2  inch.  Longitudinals  which  may  be 
I-beams,  T-beams , etc.  are  identified  by  two  dimensions  and  a 
weight.  For  example,  t.he  main  deck  longitudinals  in  Figure  3.1 
are  labeled, 

8"  X 4"  X 15  i I-T. 

The  8"  is  the  nominal  depth  of  the  section  while  4"  is  the  nominal 
flange  width.  Figure  3.2  defines  the  basic  properties  of  an  I-beam. 
The  notation  I-T  means  that  one  flange  has  been  cut  off  this  I-beam 
to  make  it  a T-beam.  Sometimes  the  I-beam  is  cut  at  the  center  of 
the  web  to  make  a shorter  T-beam.  If  this  were  done  with  the  beam 
used  in  the  example  above,  it  would  probably  be  identified  as, 

4"  X 4"  X 15  # I-T  . 

The  nominal  dimens  ons  are  then  used  to  obtain  the  dimensions  of  the 
beam  from  a handbo<-:  on  structural  design.  The  necessary  dimensions 
needed  are  the  depth  of  the  section,  the  flange  widt.h,  fiange 
thick.ness  and  web  thickness.  Another  structural  element  that  might 
be  encountered  in  reading  ship's  plans  is  the  angle,  denoted  by  the 
symbol/..  Angles  are  also  identified  by  two  dimensions  and  a 
weight.  For  example,  referring  to  Figure  3.1,  we  see  an  angle 
between  the  sheer  strake  and  the  stringer  strake  labeled 

6"  X 6"  24.2 

The  dimensions  refer  to  the  lengt.hs  of  t.he  sides  of  the  angle  and 
its  weiaht  may  be  used  to  obtain  its  cross-sectional  area  by  again 
using  the  density  of  steel.  The  area  is  equal  to  the  weight,  24.2  s. 
divided  by  the  factor  3.41.  This  will  yield  t.he  area  in  in.  sq.  . 
Thus,  in  general,  a cross-section  of  a girder  may  be  reduced  to  3 
general  elements,  i.e.,  plates,  beams  and  angles. 
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4 . APPLICATION 

4.1.  As  we  have  seen,  the  calculation  of  the  moment  of  inertia 
and  shear  area  of  a ship  cross-section,  in  principle,  involves  the 
evaluation  of  certain  integrals.  In  practice,  this  must  be  done  by 
replacing  the  integrals  by  finite  sums.  For  each  finite  element 
making  up  the  girder,  its  moment  of  inertia  about  an  arbitrary  axis, 
call  it  i,  consists  of  two  moments  of  inertia.  The  first  is  its 
moment  about  an  arbitrary  axis,  ijj,  and  the  second  is  its  moment  of 
inertia  about  its  own  axis,  iQ.  In  the  limit  of  vanishingly  small 
areas,  Iq  approaches  zero.  Consider  Figure  4.1.  Suppose  a small 
area  of  plating  is  represented  by  the  rectangular  element  of  area 
ab  rotated  by  an  angle  ^ with  respect  to  an  arbitrary  axis  3cx.  We 
can  compute  its  moment  of  inertia  about  the  rotated  x'-axis  then 
transform  this  to  obtain  the  moment  about  the  x-axis.  If  we  do 
this,  using  the  law  of  tensor  transformations,  we  find  that, 

1 ab(b^cos^0  + a^sin^iZ^)  + ab  y ^-y(ab)  (4.1) 

^NA  12  ^ 

where  y^  is  the  perpendicular  distance  from  the  arbitrary  x-axis 
to  the  centroid  of  the  area  and  y is  the  distance  between  this 
arbitrary  x-axis  and  the  true  neutral  axis.  If  we  identify, 

ij.  = ab  y^  (4.2) 

and , 

i^  = ab(b^cos^?(  + a^sin^?()  (4.3) 

we  have, 

^NA  = ^^x  ^o^  ■ (4.4) 

or, 

i^A  = i - (ab)  y ^ (4.5) 


which  is  the  Parallel  Axis  Theorem,  Equation  2.6.  The  total  moment 
of  inertia  of  the  plating  about  the  neutral  axis  is  the  sum  of  all 
the  ifj;^'s.  The  moment  of  inertia  of  beams  is  obtained  by  calculat- 
ing the  area  of  the  beam's  cross-section  from  the  four  beam 
dimensions  mentioned  previously  and  multiplying  this  by  the  square 
of  its  distance  from  the  arbitrary  xj-axis.  Its  moment  about  the 
neutral  axis  may  then  be  obtained  from  the  Parallel  Axis  Theorem. 
The  same  thing  is  done  with  angles.  Thus,  we  neglect  Iq  for  beams 
and  angles  since  they  are  quite  small.  All  of  these  calculations 
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hav’  been  reduced  to  a form  suitable  for  use  in  a computer  program. 
The  input  is  described  in  the  next  section. 

5 . PREPARATION  OF  COMPUTER  INPUT 

5.1.  As  previously  pointed  out,  the  structural  elements  that 
make  up  the  girder  may  be  grouped  into  three  categories.  These 
are  plates,  angles  and  scantlings  ( I , J , T-beams ) . The  purpose  of 
this  section  is  to  describe  how  the  data  needed  for  the  computer 
program  is  to  be  obtained.  The  tools  required  are  the  inertia 
sections,  a ruler  (preferably  divided  into  units  of  approximately 
1/50  inch),  dividers  and  a drafting  machine  for  measuring  angles. 

It  is  most  convenient  if  the  reading  operation  is  carried  out  on  a 
drafting  table  but  this  is  not  essential.  Adjust  the  machine  so 
that  the  base  line  of  the  inertia  section  (see  Figure  3.1)  is  the 
zero  degree  reference  and  the  CVK  is  the  ninety  degree  reference. 

All  angles  should  be  positive  and  between  zero  and  ninety  degrees. 

5.2.  Plating  elements  must  be  divided  up  into  rectangjes 
sufficiently  small  that  their  cross-sectional  area  is  accurately 
approximated  by  a rectangle  of  area  equal  to  length  times 
thickness.  For  decks,  bulkheads,  strakes,  etc.  that  have  little 
or  no  curvature,  a rectangle  of  much  longer  length  will  suffice 
than  for  those  which  have  a great  deal  of  curvature  and  which  must, 
therefore,  be  divided  into  rectangles  of  shorter  lengths.  The 
thickness  will  be  obtained  from  the  plating's  weight  via  Equation 
3.1.  Two  other  quantities  are  required  for  calculating  the  inertia 
of  the  plating.  These  are  the  perpendicular  distance  of  the 
centroid  of  t.he  plate  element  from  an  assumed  (arbitrary)  neutral 
a.xis  (D)*  and  the  angle  (PHI)  in  degrees  that  the  elemental  rectangle 
makes  with  this  assumed  neutral  axis  ANA,  which  is  parallel  to  the 
base  line.  Since  all  lengths  and  distances  are  in  arbitrary  units, 
scale  factors  must  be  determined  from  the  inertia  section  drawing 

in  order  to  convert  the  units  to  feet  or  inches  as  appropriate. 

This  can  be  done  by  determining  the  number  of  units  in  some  known 
distance  on  the  drawing,  for  example  the  depth  of  the  section  from 
main  deck  to  keel. 

5.3.  Angles  have  a cross-sectional  area  that  is  obtained  from 
t.he  weight  per  linear  foot  of  length  according  to  t.he  rule 
discussed  in  Section  3.1  supra.  Since  the  moment  of  i.nertia  of  the 
angle  about  its  own  neutral  a.xis  (io)  is  small  in  comparison  with 
its  moment  about  the  ANA,  only  two  quantities  are  required  for 
inertia  calculations  for  angles.  These  are  the  weight  per  linear 
foot  (WA)  (obtained  from  the  angle's  identification  on  the  drawing) 
and  the  perpendicular  distance  between  its  centroid  and  the  ANA, 
call  it  (DA).  This  distance  will  be  in  arbitrary  units.  Similarly, 
for  scantlings,  we  need  only  determine  t.he  beam's  nominal,  identify- 
ing, dimensions  and  weight  (WS) , and  the  perpendicular  distance  of 
its  centroid  from  t.he  ANA,  call  it  (DS)  . It  is  important  to  deter- 
mi.ne  the  type  (TYPE)  beam  that  is  being  considered.  These  may  be 
identified  as  I-beams,  J-beams  and  T-beams  (two  types).  .A  J-beam 

is  simply  an  I-beam  wit.h  one  end  of  t.he  flange  cut  off.  J-beams 
’'See  Table  6 . 1 for  meaning  oT  Fortran  Symbols 
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are  identified  on  drawings  by  the  symbol  J . As  previously 
mentioned,  I-beams  may  have  a flange  cut  off  to  form  a T-beam 
(referred  to  herein  as  a T- flange,  TF)  or  they  may  be  cut  in  the 
center  of  the  web  to  form  a T-beam  (referred  to  herein  as  a T-web , 
TW) . The  four  dimensions  of  the  I-beam  necessary  to  compute  its 
cross-sectional  area  are  the  section  depth  (SD) , width  of  the 
flange  (WF) , thickness  of  the  flange  (TF)  and  thickness  of  the  web 
(TV;)  . These  data  are  sufficient  for  the  calculations  carried  out 
by  the  computer  program  described  in  the  next  section. 

6.  INPUT  DATA  SPECIFICATIONS 

6.1.  The  input  data  listed  below  in  Table  6.1  is  required. 
Fortran  symbols,  formats,  units  (if  any)  and  comments  are  included. 


Table  6.1 

FORTRAN  SYr>IBOL 

FORMAT 

UNITS 

COMMENTS 

S C ALE F, SC ALE I 

2F10.0 

ft. /unit 
in . /unit 

scale  factors  for  convert- 
ing units  to  lenqth 

NP,  NA,  NS 

315 

respectively,  the  number 
of  plate  rectangles,  angles 
scantlings  that  make  up 
section 

W,L,D,PHI 

4F10.0 

lb , units 
units,  deg. 

weight/ft.  sq.  of  plate, 
length  of  elemental 
rectangles,  distance  from 
centroid  to  ANA,  angle  with 
ANA  made  by  the  rectangle 
(see  Figure  4.1) 

WA,DA 

2F10.0 

lb . , units 

weight/ft.  of  the  angle 
and  distance  from  centroid 
to  ANA 

SD,WF,TF,TW, 

DS,TYPE 

5fl0. 0 

in . , in  . , in . , 
in . , units , - 

sect,  depth,  flange  width, 
thickness,  web  thickness, 
distance  from  centroid  to 
ANA,  beam  type 

6.2.  A few 

additional 

comments  are 

in  order.  First,  the 

determination  of  the  centroid  of  the  plate  rectangles,  angles  and 
scantlings  is  not  critical.  On  some  inertia  sections,  the  centroid 
of  scantlings  is  indicated.  If  it  is  not,  its  position  may  be 
estimated  by  considering  how  it  will  be  shifted  by  the  removal  of 
either  a flange  or  part  of  the  web  and  the  flange.  Second,  when 
reading  inertia  sections,  it  is  useful  to  prepare  data  sheets 
appropriate  for  recording  the  requisite  data.  These  should  be 
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labeled  with  the  section  number,  ship,  date  and  number  of  plate 
elements,  angles  or  scantlings.  Finally,  it  should  be  noted  that 
although  the  computer  program  is  set  up  to  read  data  for  an  entire 
ship  cross-section,  it  is  necessary,  because  of  the  symmetry,  to 
read  only  one-half  of  the  cross-section.  When  digitizing  the  data, 
it  is  necessary  to  pxinch  one  card  and  then  this  card  may  be  dupli- 
cated. Remember,  however,  that  in  a few  instances  there  may  be 
P and/or  S asymmetries  in  the  plating.  This  has  been  discussed  in 
Section  3.2. 

7.  COMPUTER  PROGRAM  AND  OUTPUT 

7.1.  The  computer  program  that  carries  out  the  longitudinal 
strength  calculations  is  written  in  Fortran.  A listing  is  provided 
as  Appendix  A.  All  necessary  scaling  is  done  internally  so  that  it 
is  unnecessary  to  convert  lengths  and  distances  to  real  lengths  and 
distances  externally. 

7.2.  The  program  computes  the  moment  of  inertia  and  static 
moment  about  an  arbitrary  assumed  neutral  axis  and  then  uses  the 
Parallel  Axis  Theorem  to  obtain  the  moment  about  the  true  neutral 
axis.  It  also  computes  the  total  cross-sectional  area  of  the 
strength  elements  that  make  up  the  section.  The  calculations  are 
carried  out  in  stages.  First  the  contribution  of  plating  is  consi- 
dered, then  the  angles  are  added,  then  the  effect  of  the  scantlings 
is  considered.  The  cross-sectional  area,  moment  of  inertia  and 
static  moment  are  printed  out  at  these  three  stages.  This  is  useful 
in  compa'''ing  the  relative  contributions  of  these  types  of  structural 
elements.  Also  printed  out  is  the  distance  between  the  true  neutral 
axis  and  the  assumed  neutral  axis.  The  cross-sectional  area,  moment 
of  inertia,  static  moment  and  distance  between  true  and  assumed 
neutral  axis  are  denoted  respectively  by  SUMA,  SUMI , SUMQ,  and  DG. 
After  obtaining  the  final  totals  for  these,  the  program  computes  the 
thickness  of  plating  on  the  neutral  axis  (H) , the  moment  of  the 
area  above  the  neutral  axis  (MA)  and  uses  Equation  2.10  to  obtain 
the  shear  area  (ASHR) . The  units  of  output  quantities  are  given 

in  Table  7.1. 
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OUTPUT  VARIABLE 
L,T,A,D,PHI, 10 


AA,DA 

AS, DS, TYPE 
SUMA,SUMQ,SUMI 

DG, INA 

SA,SQ,MA 

H,ASHR 

A typical  data 
program  listed 


Table  7.1 


DIMENSION 

in . , in . , in . sq . , 

ft . , deg . , in . sq . - f t . sq . 


in . sq. , ft . 


in.  sq. , ft. 


in.sq. ,in.sq.-ft.  , 
in.sq.-ft.sq. 


ft. , in. sq. -ft. sq. 


in.sq. , in . sq. -ft . , 
in. sq.-ft 


in . , in . sq . 


set  is  included  in  Appendix 
there . 


COMMENTS 

The  length,  thickness 
and  area  of  the  elemen- 
tal rectangles  that 
make  up  the  plating, 
distance  of  centroid 
from  ANA,  the  angle 
made  by  the  rectangles 
with  respect  to  the 
ANA,  and  the  moment  of 
inertia  of  the  rectangle 
about  its  own  axis . 

cross-sectional  area 
of  the  angle  and 
distance  of  centroid 
from  ANA 

cross-sectional  area 
of  the  scantling, 
distance  of  centroid 
from  ANA,  type  of  beam 

cross-sectional  area 
of  the  structural 
elements,  static  moment 
about  ANA,  moment  of 
inertia  about  ANA 

distance  between  ANA 
and  true  neutral  axis, 
moment  of  inertia  of 
section  about  the  true 
neutral  axis 

area  and  static  moment 
of  area  above  ANA, 
moment  of  area  above 
true  neutral  axis 

thickness  of  plating 
on  the  neutral  axis, 
shear  area  of  cross- 
section 

for  use  in  the  computer 
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FIG.  2.1  PLANE  AREA  WITH  CENTROID  C 
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FLANGE  THICKNESS 


SECTION  DEPTH 


FIG.  3.2  NOMINAL  BEAM  DIMENSIONS 
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APPENDIX  A 


or MF NS I ON  A(200).L(200)»0<POO),lO(?00),r(POO)«PHl(200)»W(200) 
or  MENS TON  wA(200).DA(200).TYPE(20  0).AA(20  0)tASC(200) 
dimension  so (200). wF (200). TF(200)*TW(200) *05(200) 

HEAL  L.IO.INA.MA 
data  PT/l.lAlSq?f./ 

REAP  (‘^.10)  SCALFF, SCALE! 

REAP (5. 10 ) YANA 
yana=yana*scalef 

C YANA  IS  the  OISTAmcE  FROM  THE  KEEL  TO  THE  ASSUMED  NEUTRAL  AXIS 

C SCALE  CONVERTS  ARPITRARY  UNITS  INTO  FEET  OR  INCHES 

C A IS  THE  AORAY  OF  ELEMENTAL  AREAS 

C L IS  THE  ARRAY  OF  ELEMENTAL  LENGTHS 

C T IS  the  array  OF  elemental  THICKNESSES 

C 0 IS  THE  ARRAY  OF  DISTANCES  BETWFEnj  The  ELEMENTAL  CENTROIDS  AND  Th£ 

C ASSUMED  neutral  axis 

C lO  IS  THE  ARRAY  OF  MOMENTS  OF  THE  ELEMENTAL  AREAS  ABOUT  THEIR  OvKN 

C AXIS  parallel  TO  THE  ASSUMED  NEUTRAL  AXIS 

C Phi  IS  THF  ARRAY  OF  ANGLES  THAT  EACH  ELEMfcNTAL  LENGTH  MAKES  wITh  ThE 

C assumed  NFUTRAL  AXIS 

C w IS  THF  array  of  plate  'WEIGHTS  ASSOCIATED  WITH  EAC”  ELEMENTAL  AREA 

READ  (S.20)  NP.NA.MS 

NP  IS  THE  NUMPER  OE  PLATE  ELEMENTS  TO  PE  PEAQ  IN  AS  OATA 

NA  IS  THE  NUMBER  OF  ANGLES  TO  8E  READ  IN  AS  OATA 

NS  IS  the  NUMHER  of  SCANTLINGS  TO  PE  READ  IN  AS  DATA 

REAO  (S.IO)  (w ( I ) ,L ( I ) ,D ( I ) .PHI ( I ) , T=1 ,NP) 

SUMA=0 . 

SUM(5  = 0 . 

SUMlrO. 

SA  IS  the  SUM  OF  THE  AREAS  LYING  ON  OR  ABOVE  THE  ASSUMED  NA 
SO  IS  THE  SUM  OF  STATIC  MOMENTS  OF  ARFAS  LYING  ON  OR  ABOVE  ThE  NA 
SA  = 0. 

S0  = 0. 

WRITE  (B.2S) 

WRITE  (6.211 
WRITE (6. 2S) 

DO  40  1=1. MP 
T ( I ) =W ( I ) /40.H 
L(I)=L(I)*SCALEI 
A ( I ) =L ( I ) *T ( I ) 

D(I)=0(I)*SCALEF 
SUmA=SUMA*A ( I ) 

SUMO=SI)MO*A  ( I ) *0  ( I ) 

TH=PHI (I)*PI/180. 

10  ( I ) =A ( I ) /I  728.*  C (T ( I ) *COS (TH)  ) ♦ ( L ( I ) *S I N ( TH ) ) **2) 

SUMI=SUMI*A ( I ) *D ( I ) *0 ( I ) *10 ( I ) 

WRITE  (6. SO)  L(I)  .T(I) .A(l) .0(1)  .PHI (1) ,10(1) 

40  continue 

WRITE  (6,27) 

WRITE (6.277) 

WRITE  (6,50)  SUMA.SUMO.SUMI 
DG=SUMO/Sl'MA 
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r Dh  IS  the  distance  between  the  assumed  neutral  axis  and  The  true 
c neutral  axis 

INA»SUhI-SUMA*OG*DG 

C IT  IS  ASSUMED  that  the  SHIP  IS  SYMMETRIC  ABOUT  ITS  CENTER  LINE 

■PITF  (6,?B) 

write 

write  (A, SO)  DG»INA 

IE  inA.EO.O)  go  to  7S 

READ  (S.IO)  (WA  ( I ) ,DA ( I ) , Isl «NA) 

wA  IS  the  array  oe  weights/ft  for  the  angles 

DA  IS  the  array  of  distances  between  the  angle  centroids  and  the 

assumed  neutral  axis 

WRITF  (A,?4) 

write  (A. 37) 

WRITE (A, 377) 

DO  70  I*1»NA 
AA(I)sWA(I)/3.410R 
AA  IS  THE  ANGLE  AREA  IN  IN  SO 
DA ( I ) =DA ( I ) *SCALEF 
SUMAsSUMA*AA  ( I ) 

SUMOsSUMG* AA ( I ) «DA ( I ) 

SUMI=SUMI*AA ( I )*DA ( I ) *DA ( I ) 

WRITE  (A, SO)  AA{I),DA(I) 

70  continue 

WRITE  (A.?7) 
write (A»277) 

WRITE  (6, SO)  SUMA.SUMO.SUMI 
DG=SUM0/SUMA 

ina=sumi-suma*dg*dg 
write  (6,?B) 

WRITE (A. 288) 
write  (A, 50)  DGtiNA 

7S  READ(5.12)  (SD ( I ) .WF ( I > ,TF ( I ) »TW { I ) ,DS fl ) »TYPE ( I ) . 1=1 »NS) 

SD  IS  the  array  of  scantling  depths 
WF  IS  the  array  of  flange  widths 
tf  is  the  array  of  flange  thicknesses 
Tw  IS  THE  array  of  WEB  THICKNESSES 

DS  IS  THE  ARRAY  OF  DISTANCES  BETWEEN  SCANTLING  CENTROIDS  AND  THE 
ASSUMED  neutral  AXIS 

type  identifies  The  scantling  as  i.j.  web  cut  t.  or  flange  cut  t 

WRITE  (A.2A) 

WRITE  (6.38) 

WRITE (A. 388) 

DO  100  1=1. NS 
DS(I)=DS(I)*SCALEF 


IF 

(TYPE ( I ) 

.EQ.2HI  ) 

GO 

TO 

105 

IF 

(TYPE(I) 

.E0.2HJ  ) 

GO 

TO 

no 

IF 

(TYPE  (I  ) 

.E0.2HTW 

) GO 

TO 

115 

IF 

(TYPE ( I ) 

.EQ.2HTF) 

GO 

TO 

120 

IL 
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1 OS  AS  = SD ( I ) »WF  < I ! - (SD ( I ) -2,*TF ( I ) ) * (WF ( I ) -TW ( I ) ) 

GO  TO  130 

110  AS  = SD(I)*WF<I)-.5*(WF(I)-TW(I) ) * ( 2 , *50 ( I ) -3 . *TF  ( I ) ) 

GO  TO  130 

1 1«^  AS=.5*  <SO(  I ) *WF  (I ) - (SO  ( I ) -2,*TF  (I)  ) * (WF  ( I ) -TW  ( I ) ) ) 

GO  TO  130 

1 20  AS  = SO ( I ) *WF ( I ) - (SD ( I ) -2.*TF ( I ) ) * (WF ( I ) -TW ( I ) ) -WF ( I ) *TF ( I ) 

130  ASC(I)=AS 

SUMA  = SUMA*ASC ( I ) 

SUM0  = SUM0*4SC ( I ) *DS  ( I ) 

SUMirSUMI  + ASC ( I ) *05 ( I ) *05  ( I ) 

WWITF  (6,SS)  AS*OS ( I ) tTYPE  ( I ) 

100  CONTINUE 

WRITE  (6.27) 

WRITE (6,277) 

WRITE  (6,50)  SUMA,SUMQ,SUMI 

DG=SUMO/SUMA 

INA  = SUN'I-SUMA*06*0G 

WRITE  (6,20) 

WRITE (6.288) 

WRITE  (6,50)  DG.INA 

CALCULATION  OF  THICKNESS  OF  PLATING  ON  THE  NEUTRAL  AXIS  AND  THE 
Shear  area  of  the  section 

H=0. 

DO  36  1=1, MP 
TH=PHI ( I ) *PI/180. 

DEL=.5*L(I)*SIN(TH)/12. 

IF  (D(I)-OG  ,LT.  0.  .and.  ABS(D(I)-0G)  .GT.  DEL)  GO  TO  36 
IF  (A8S(0(I)-OG)  .GT.  DEL)  GO  TO  35 
H=H+T ( I ) 

A ( I ) =T  ( I ) *( .5*L( I ) ♦ (D( I ) -OG) *12./SIN(TH)  ) 
D(I)=.25*L(I)*SIN(TH)/12.*.5*(D(I)-0G) 

35  SA=SA^A(I) 

SO=SOfA ( I ) *0 ( I ) 

36  CONTINUE 

IF(H  .NE.  0.)  GO  TO  39 

NP1=NP-1 

K=1 

DO  1 1=1, NPl 

IF (APS (PHI ( I ) ) .LT.  10.)  GO  TO  1 
0=AP5 (0(1) -DG) 

IF(0  .LE.  APS (0 (K) -DG) ) K=I 
1 CONTINUE 
H=T (K) *2. 

A(K)=T(K)*(.5*L(K)*(D(K) -OG) *12./SIN (TH) ) 

0 (W) =,PS*L (K) *SIN(TH) /12.*.5* (D (K) -OG) 

SA=SA*?.*A  (K ) 

SQ=SQ+?.*A (K) *0(K) 

39  continue 
39  continue 
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IF  (NA.FQ.O.)  60  TO  133 
DO  42  1 = 1, 

IF  (HA (I ) -06.lt. 0. ) GO  TO  42 
SA  = SA* AA  ( I ) 

Sij=SO*  AA  ( I ) *OA  ( I ) 

42  CONTInUF 
133  00  F7  1=1, MS 

IF  (ns ( I ) -06.lt. 0. ) GO  TO  57 
SArSA* A5C ( T ) 

S(0  = SQ  + ASC  ( I ) *DS  ( I ) 

^^7  CONTINUE 

MA  = so-no*';A 
ashr  = 2.*ina*h/'-ia 

r A5hR  IS  THE  SHEAR  AREA 

AShR=ASHR«1 2. 

c ashr  converted  to  SQ.IN. 

D6=YANA  ♦ OG 

»'RITF  (F,?2)  DG.H,SUMA,ASHR,  INA 

in  format  (4F10.0) 

12  f-ORMAT  (5F10.0,A2) 

20  format  (315) 

21  format  (//5X«LENGTH*7X*THICKNESS*R)(*AHEA*fcX*0lST.  FROM  NA*6X*ANGlE* 
17X»I0  OWN  AXIS*) 

2?  format  (/IX,  2RHDISTANCE  FROM  KFEL  TO  N.A.  = »FS.2,  5H  FEET/lX.pRH 
IThICKNESS  of  PLATE  ON  N.A.  = »F8.3,7M  I NCHE S/ 1 X . 29H ARE A OF  CROSS  S 
lECTlON  = F5.0.11H  SO.  I NCHE S/ 1 X » 3 1 HSECT I ON  SHEAR  AREA 

1 = , FH.O.llH  SO.  IMCHES/1 X.29HSECTI0N  MOMENT  OF  INERTIA  = 

1 . FR.O.IRHSO.  FEET-SO.  INCHES) 

23  format (//60X*PLATING*//) 

24  format  {//f.nx«AN6Lf  S*//) 

?5  F0RMaT(4X«(ImChES)*7X*( INCHES) *7 X* (SO.IN.)«7X*(FEET)*8X* (DEGREES) « 
13X*(In.S0.-FT.SQ.)*/) 

26  format (//60X*SCANTLINGS*//) 

2 7 format (//5X*AREA*1 0 X*M0MENT*RX* INERT  I A* ) 

277  format (3X*(SQ.IN.)*5X*(IN.S0.-FT.)*2X*(IN,S0.-FT.S0.)*/) 

2«  format (//I X*niST.  FROM  NA*3X*I  about  NA*) 

2M9  F0RMAT(6X*(FT)*5X*(IN.SQ.-FT.S(3.)*/) 

17  format (//5X*AREA*7X*DIST. FROM  NA*) 

37  7 format  (3X*(TN.SC).)*10X*(FT,)*/) 

IS  format (//SX*AREA*7X*0TST, FROM  NA*8X*TYPE*) 

38  B format  (3X*(S(3.IN.)*RX*(FT.)*/) 

so  format (Fin,?,5X.F11.3.4X,F11.3f4X»F10,2»5X»F10.2»5X.F12.4,3X) 

55  format (FI 0.3.5X »F ) 1 . 3. 1 1 X , A2) 

STOP 

END 
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0. 
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0. 
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0. 

2“. 

9.6 
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0. 

?9. 

9.6 
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0. 
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13. 
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0. 
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9.6 

77. 

1 1 . 

-130. 

9.6 

27. 

10. 
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16. 
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16. 
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